arXiv:astro-ph/0512641v2 20 Jan 2006 


Analytic Results for a Flat Universe Dominated by Dust and Dark Energy 
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We find the solution for the scale factor in a flat Universe driven by dust plus a component 
characterized by a constant parameter of state which dominates in the asymptotic future. We also 
present an analytic formula (in terms of hypergeometric functions) for the past light cone in such 
a universe. As applications of this result, we give analytic expressions for the Luminosity Distance 
and for the Acoustic Scale, where the latter determines the peaks positions in the pattern of CMB 
anisotropies. 

PACS numbers: 98.80.-k,98.80.Es 


I. INTRODUCTION 

The past light cone, defined as the comoving distance 
r travelled by a light signal from time t to the present 
time to (> £), in cosmology is given by JjL] 

r(t ' to)= iF7m s, ‘( HaVWS l $>)’ (1) 

where a(t) is the scale factor, H 0 is the Hubble pa¬ 
rameter (defined as H = din a/dt) computed at to and 
= — K/{HqOq) is the present curvature energy den¬ 
sity (K = 1,0,—1 for closed, flat, open spatial sections) 
with ao = a(to). Depending on K = 1,0,—1 respec¬ 
tively, the function Sk- 

{ sin(x) 
x 

sinh(a;) 

contains information on the geometry of space. 

For flat spatial sections Eq. 0 simplifies to 

r(t, t 0 ) = J dt'/a(t'). (2) 

Changing variable of integration from t to the redshift z 
(defined as z = 1/a — 1, with a(t o) = 1), one obtains 

r ( z ) = f dz'/H(z'). (3) 

Jo 

Since H(z) is given by the Friedmann equation (see 
Sec | 0 , Eq. 0 allows the computation of r without solv¬ 
ing for a(t). Of course, once r(z) is known, one still has 
to solve the Friedmann equation in order to obtain the 
dependence on time. 

The integrations 0.0 are the starting points for the 
computation of some distances of astrophysical interest 
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as the Proper Distance (or Horizon), the Angular Diam¬ 
eter Distance and the Luminosity Distance j2j. The co¬ 
moving distance 0 is also needed for the computation 
of the Acoustic Scale that is the characteristic angular 
scale of the peaks in the Power Spectrum of the Cosmic 
Microwave Background (henceforth CMB) anisotropies. 

An analytic form for r is known for few simple cases 
lj because is difficult to perform the integration either 
in form of Eq. 0 or Eq. 0 . In the case of dust plus 
curvature, the time evolution for a and r are known 0. 

Here we present the analytic solution of the scale fac¬ 
tor and the exact computation of the past light cone in 
a flat Universe driven by Dust (whose energy density 
pc(t) oc a -3 ) and Dark Energy ( px(t ) oc a a ) with a 
parameter of state (i.e. pressure over energy) constant in 
time. The important case of ACDM is obtained setting 
a = 0. We apply this formula to the following cases: first, 
we provide the exact expression of the Luminosity Dis¬ 
tance (and we give an expansion of the past light cone at 
large z) and then we write down the analytic expression 
for the Acoustic Scale, relevant for the peaks analysis of 
the angular power spectrum of CMB anisotropies. 

The paper is organized as follows: in section[0we solve 
the Friedmann equation introducing a new time coordi¬ 
nate. In section EH we find the analytic formula for the 
past light cone. We apply our results to the luminosity 
distance in section HYI and to the acoustic scale in section 
El We conclude in section eh 

II. SOLUTION OF THE FRIEDMANN 
EQUATION 

In a flat universe, considering the FRW metric ds 2 = 
—dt 2 + a 2 (t)dx 2 , filled with dust and another component 
with a constant parameter of state (w x = —1 — a/3) we 
have 

pc(t) = pca~ 3 (t) and p x (t) = p x a a (t) , (4) 

where p x , pc do not depend on time and are the en¬ 
ergy density of the two components at present time. 
For a > —2 the A—component drives the universe into 
asymptotic acceleration and for a > 0 such component 
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violate the dominant energy condition (wx < — 1). The 
scale factor a(t) satisfies the Friedmann equation 


H 2 



8t tG 


Pxa 



(5) 


The total state parameter is 


wtot(o) 



p x a 


£*+3 


Px a 1 


ck+3 . 


PC 


( 6 ) 


and interpolates between 0 and wx ■ It is not easy to find 
solutions for a two-component system (see however the 
radiation-matter universe and one of its generalizations 
0), but Eq. © can be integrated in the a = 0 case (i.e. 
ACDM model) When a / 0, it is no more possible 
to obtain analytically a = a(t, a ) but only its inverse 
t = t(a,a). This problem can be solved defining a new 
time coordinate r: 


dt = dr a “/ 2 , 


(7) 


such that in this new coordinate, eq. © is given by 


Hi = 


87 tG 


Px 


Pc 


7 3-t-c* 


( 8 ) 


with H t = a'/a, where the prime stands for derivative 
with respect to the time r. Note that only for a = —2, i.e. 
when the X-component behaves like a curvature term, r 
becomes the conformal time 77 (defined as dp = dt/a). 
Eq. © can be integrated for a > —3 with the result: 


a(r) = 
x sinh 


ru 


1 — f lc 
3 + oc 


1/(3+a) 


(1 - n c ) 1/2 H 0 T 


1 2/(3-f-a) 


(9) 


where f lc = 8 ttGpc/3H 2 = 1 — f lx is the present dust 
density ratio 0. For r —> 0, we find the standard dust 
behaviour a ~ r 2 ^ 3+a ^ ~ t 2 / 3 . Notice that when a —> 0, 
Eq. © recovers the known solution for ACDM model 0. 


III. PAST LIGHT CONE 

In order to perform the integration ©, we make use 
again of Eq. © so that is possible to rewrite the comov¬ 
ing distance as 


where 1/A = (— 1 + fl c 1 ) 1 / | ' 3+Q \ The integration left in 
Eq. ITT1) can be performed giving (see 15.3.1 of |6|) 


r{z) 


1 2 
H~on^ 



1 + z) ’ fl c (l + Z ) 3+a J 


,( 12 ) 


where the definition of the redshift z has been used and 
where / is the following hypergeometric function (6] 




2Fl |_2(a + 3)’ 2’ 1 + 2(a + 3) ’ 

(_ y )ITO / 1 1\ 

(6 + 2a) v \6 + 2a 2J ’ 


(13) 


and in the last equality we have also rewritten the hyper¬ 
geometric function in terms of incomplete Beta function 

B x { 7,/?) li¬ 
lt is interesting to expand the past light cone for large 
z, obtaining 


-(*) = 




1 — f lc 


n 


c 


2 Hi 

1 


1 


(1 + z) 1 / 2 


2(a + 3) ’ 2 

+ 0 


I.i + 


2(a + 3) ’ 


(1 + z)“+ 7 / 2 


, (14) 


where the leading term (that is constant with respect 
to the redshift z) is the conformal age of the universe 
already given in :7|. Note that the ACDM cosmic age in 
cosmic time is given in lj. 


IV. LUMINOSITY DISTANCE 

Eq. m is the starting point for the computation of 
distances of astrophysical interest such as the Luminosity 
and the Angular Distances. The Luminosity Distance, 
dL, and the Angular Distance, in a flat universe are 
given by 


d L (z) = (1 + z) r(z ), 

(15) 

d A {z) = (1 + zy 1 r(z ). 

(16) 


r(a) 


da 


a l+a/2 a ' ‘ 


( 10 ) 


By making explicit the form of H T in the integral, 
Eq. (ITTift reads 


Let us focus on the Luminosity Distance. From 
Eq. (1151) . dL can be expanded for small z by using its 
form in terms of incomplete Beta function 

HodL(z) = z + d 2 Z 2 + c? 3 Z 3 + djz 4 + 0{z 3 ), (17) 


"(a) = 


(c + 2) 

1 - n c \ 

flc 

■1/A 


H 0 (l - Dc) 1 / 2 


dx 


• a/A 


c 1 / 2 (l 


>+3+0: 


)l/2 


( 11 ) 


with 


d2 


(4 + a) - (3 + a)Q.c 
4 


(18) 
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FIG. 1: Left Panel: Exact luminosity distance Hq dj, (solid line) and approximated luminosity distance Hod app (up to third 
order, dotted line), both vs redshift 2 . Right Panel: p vs z where p is implicitly defined by dn = (1 + p)d apP where dr, is the 
exact luminosity distance and d app is the approximated (up to third order) luminosity distance. In both panels Qc = 0.3 and 
a = 0 . 


d 3 


^ [a(4 + a) - 2(3 + a)(5 + 2a)fl c + 

+3(3 + o:) 2 f2^] , (19) 


+3(3 + a) 2 (20 + 9a)flc - 15(3 + a) 3 !^] (20) 


d,4 


yyy [(-2 + a)a(4 + a) 

— (3 + a) (40 + a(50 + 13a))O c 


Once a is set to zero, d ,2 and d 3 recover the known co¬ 
efficients in fiat ACDM model |8) with the deceleration 
parameter qo = £lc/2 — fix = 3flc/2 — 1 and with the 
jerk parameter given by jo = f lc + fix = 1- 


J 


In Fig. [I] we show the comparison between the exact 
(see Eq. (1151 1 and the approximated (up to the third 
order, see Eq. 03 ) Luminosity Distance as a function of 
z for fl c = 0.3 and a = 0. For istance at z = 1.5, an 
error of 9.8% is made, if one considers the approximated 
instead of the exact luminosity distance. 

In Fig. [3 we compare the exact luminosity distance 
with the first three approximated luminosity distances in 
a flat ACDM model. We note that the approximation 
of order z 2 always overestimate the exact luminosity dis¬ 
tance while the approximations of order z 3 and z 4 always 
underestimate the exact luminosity distance. From this 
analysis it is clear that already for z ~ 2 the use of a third 
order expansion may flaw the comparison of observations 
with theoretical predictions. 

In the left panel of Fig. 0 we set flc = 0.3 (and hence 
fix =0.7) and plot Eq. (tTSll for different values of a. We 
notice that the dependence on a is quite weak. In the 
right panel of Fig. 0 we set a = 0 (i.e. ACDM model) 
and plot Eq. 03 for different values of the ratio fix/f 1c- 
Of course larger is this ratio and larger is the luminosity 
distance for each fixed redshift. We notice that the de¬ 
pendence on this ratio is stronger than the a dependence. 


V. ACOUSTIC SCALE 


The characteristic angular scale of the peaks of the 
ang ular power spectrum in CMB anisotropies is defined 
as 


_ ‘I's^Zdec) 

^A ^(^-dec) 


( 21 ) 


where r s (zd e c ) is the comoving size of the sound horizon 
at decoupling, r(zd e c) is the comoving distance at decou¬ 
pling given in Eq. m and £a is the multipole associated 
to the angular scale 9a, also called acoustic scale. In fact, 
the expression d misses the contribution coming from 
the radiation (fl ra d)- A numerical study shows that we 
get an overestimate for £a that depends on a and in the 
worst case it is of the order of 0.5%. The comoving size 
of the sound horizon is given by 

~ z dec) An 

r s (zdec) = / - TT —c s (a), (22) 

J Q a z H{a) 

where c s is the sound speed given by 0 


c s (a) = 1 /^3+ (23) 

where ST^ and fl 7 are respectively, the present density 
ratio for baryons and for photons. The integral d can 
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FIG. 2: Exact luminosity distance (solid line) and three ap¬ 
proximated luminosity distance for flc =0.3 and a = 0 vs 
redshift. Long, short dashed and dotted line are second, third 
and fourth order in z expansion, respectively. Future cosmo¬ 
logical interpretations of high-redshift data )9] need the accu¬ 
racy of the exact solution. 


be analytically performed once fix , that enters in H (a), 
is set to zero. The result is m 


T s (& dec) 


4 / Qrad 

3-ffo V flc^lb 

, \/l 1 Rdec T \J Rdec T Req 

X m - - - - 

1 + \J Req 


,(24) 


with R(a) = 3(flb/(4n 7 ))a and where the label dec 
stands for “computed at decoupling” while eq stands for 
“computed at equivalence” (between radiation and mat¬ 
ter). A numerical study shows that in the worst case this 
approximation (i.e. fix = 0 in r s ) gives an error of the 
order of 10 -5 %, then completely negligible. 


We are now ready to make a prediction for the acoustic 
scale £ a defined in Eq. <EU). We consider a = 0 (i.e. 
ACDM) and set a dec = 1/1089, fl c = 0.27, fl b = 0.046 
and fix = 0.73. In order to compute z eq we made use of 
the following relation 0 

5464(^/0.135) 

69 (T C mb/2.725) 4 (1 + p,/p 7 ) ’ [ } 

where cu m = cu b + w c j m = flch * 2 (with h = 
ido/(100kms _1 Mpc -1 ), p v /p 1 is the ratio of neutrinos 
to photons densities and Tqmb is the CMB temperature 
in °K. Setting Tcmb — 2.725, h = 0.72, p v lp 1 = 0.6851 
E3 we get z eq = 3362. Once z eq is known, it is possible 
to compute fl 7 through 


fi 7 = 


flc 

(1 + Pvj P^){ z eq + 1 ) 


(26) 


obtaining fl 7 = 4.76 10 -5 . Thus we have all we need 
to compute the acoustic scale defined in Eq. (12 111 . We 
obtain 


i A = 301.76, (27) 

which is in agreement with the Wilkinson Microwave 
Anisotropy Probe (WMAP) lcr result [l(J]: (WMA.P _ 
300 ± 3. 

So our analytic calculation of the acoustic scale agrees 
with the value reported by the WMAP team. This is 
a further application of our results to a different set of 
observations, i.e. CMB anisotropies, which require a cor¬ 
rect evaluation of the past light cone at high redshift, 
where a perturbative expansion does not work and only 
numerical methods were used so far. 

Therefore, by using Eq. m and Eq. o, the acoustic 
scale is given by: 


J 


37 t fa 2 F, [- 1 /( 60 ^), 1/2,1 - 1 /( 60 ;*), - (1 -fl c )l^c] - 1/(1 + Az e c) 1/2 

2 y flj In [i/l + Rdec + \JRdec + Req\ — In [l + ^ R e q\ 


where the dependence on the cosmological parameters 
and in particular on lux is explicit. Our Eq. 12811 is com¬ 
pletely analytic compared to A£ given in [l2(, in which 
r{zdec) contains a function F that is computed numer¬ 
ically 0 and r s (z dec ) contains an averaged speed of 
sound that is fixed to a particular numerical value. In 
Ref. 0 the dark energy parameter of state is a function 
of time, but is taken as a time-average for the final expres¬ 
sion of A£ in the models considered in 0- In Fig. 0 we 
show the dependence of the acoustic scale on a, keeping 
the other parameters fixed to the values reported above 
in the text. 


VI. CONCLUSION 

We have found the analytic solution for the scale factor 
and past light cone for a flat universe driven by dust and 
another X component with a parameter of state constant 
in time which dominates in the future. If wx < —1/3, X 
is the simplest possible parametrization for Dark Energy. 
While Eq. (0 generalizes the previously known solution 
for ACDM || to a / 0, the result (t*H21) extends Q for 
finite z, which is important for accurate quantitative pre¬ 
dictions. 
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FIG. 3: Left Panel: Exact luminosity distance vs redshift for Qc = 0.3 and a = —2 (solid line), a = —1 (dotted line), a = 0 
(short dashed), a = 1 (dashed line) and a = 2 (long dashed line). Right Panel: Exact luminosity distance vs redshift for a = 0 
and for Qx/^c = 9 (solid line), Tlx/Tic = 7/3 (dotted line), Tlx/Tic = 1 (short dashed line), f lx/Tic =3/7 (dashed line), 
fix /Tic = 1/9 (long dashed line). 
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FIG. 4: Acoustic Scale vs a. Our analytic evaluation of 
the acoustic scale lies between the dashed lines. The higher 
dashed curve is computed with Eq. m while the lower is the 
99.5% of the previous one. In this way we give an estimate 
of the order of magnitude of the error it is made, neglecting 
the contribution of the radiation in Eq. Q . Horizontal lines 
represent the lcr contour to the value given by M- 


Eq. U3 has been applied to the following cases: the 
computation of the luminosity distance, relevant for su¬ 
pernovae, and the acoustic scale, which determines the 
position of the peaks in the pattern of CMB anisotropies. 
Both of these results are original. For the second case, 
our analytic prediction is in agreement with the value 
reported by the WMAP team. Note that this latter ap¬ 
plication requires a correct evaluation of the past light 
cone at high redshift, where a perturbative expansion 
does not work and only numerical methods were used so 
far. 

A parameter of state constant in time cover a wide 
range of cosmological models which include dust. From 
a fundamental physics point of view, X corresponds to 
a cosmic string network or a curvature term when a = 
—2 or to a domain walls background when a = — 1. A 
quintessence model with an exponential potential which 
does not have tracking behaviour leads to a component 
with — 2 < a < 0 (for a > 0 one has to reverse sign of 
the kinetic term of a canonical scalar field, i.e. consider 
a phantom with an exponential potential). 
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In analogy this trick can be used to find the solution for a bouncing universe by replacing px with —px, generalizing the 
matter plus positive curvature solution 0. In terms of the new time r the solution for the scale factor is: 


a(r) 



i/(3+°0 


sin 


f 3 + 0! 

I 2 



r + ip 


2/(3+a) 


where p is taken 0 for a > —3 and n/2 for a < —3. In this way the bounce is always at r = 0. Note that the allowed values 
for a are — 3 + 1/n with n being an non-vanishing integer. 
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